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Abstract 

| The CP asymmetry in the B — ► KsK$Ks decay is being measured by the two B factories. A 

large deviation of the CP asymmetry Sk s k s k s horn —S^k s and/or of Ck s k s k s from zero would 

> ' 

imply new physics in b — > s transitions. We try to put upper bounds on the Standard Model size 
of these deviations, using SU(3) flavor relations and experimental data on the branching ratios 
of various decay modes that proceed via b —* d transitions. We point out several subtleties that 
distinguish the case of three body final states from two body ones. We present several simple 
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Qh! relations that can become useful once all relevant modes are measured accurately enough. 
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I. INTRODUCTION 

The Belle and Babar [2] experiments have recently presented their first results on the 
CP asymmetries in B — ► KsK$Ks decays. The average of the two measurements is given 
by 3 

Sk s k s k s = -0.26 ±0.34, 

(1) 

Ck s k s k s = -0.41 ± 0.22. 
The B — > KsKgKs decay is a flavor changing neutral current process and, consequently, 
does not proceed via tree level diagrams. Within the Standard Model, the b — > s penguin 
contributions are dominated by a single weak phase, that is the phase of V* b V cs . The effects 
of a second phase, that is the phase of V* b V us , are CKM suppressed by 0(\ 2 ). Neglecting 
he latter contributions, and taking into account that the KgKsKg state is purely CP even 



4] 



the SM predictions are then as follows: 



Sk s K s K s ~ S^Ks 

C Ks k s k s ~ 0, (2) 

where, experimentally, S$k s — 0.726 ± 0.037 |3j. These predictions are valid also in exten- 
sions of the SM where the B° — mixing amplitude is possibly affected by new phases, 
but the b — > s decay amplitudes are not. 

In order to understand whether violations of eq. (J2J) signal new physics, it is necessary to 
estimate or, at least, put an upper bound on the CKM-suppressed SM contributions. Such a 
calculation involves, however, in addition to the CKM factors, hadronic physics. Currently, 
no first principle method for calculating hadronic matrix elements has beenproven to work 



have 



to a high level of precision. Furthermore, existing methods (for example, 
only been applied to two body final states, while our interest here lies in the three body 
mode B — > KsKsKs- 

In this work we use the approximate SU(3) of the strong interactions to constrain the rel- 
evant hadronic matrix elements i9] . While this method has the advantage of being hadronic- 
model independent, it has the following two weaknesses. First, SU(3) breaking effects can 
be of order 30%, so that our results cannot be trusted to better accuracy than that. Sec- 
ond, since we have no information about the strong phases that are involved, we make the 
most conservative assumption, whereby all amplitudes interfere constructively. This leads 
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to upper bounds on the deviations from eq. (J2J) that are often much weaker than the ac- 
tual deviations expected in the Standard Model. Furthermore, the quality of our upper 
bounds depends on the precision of current measurements. Thus, one should not think of 
our bounds as estimates of the deviation expected within the Standard Model. They are 
only approximate (to 0(0.3)) and, in most cases, very conservative upper bounds (with the 
advantage of being model independent). 

This type of analysis has been previously applied to CP asymmetries in decays into two 
body final states u sing the full SU (3) symmetry or into three body final states using an 



SU(2) subgroup [12, LL3J, |1J] . For the mode of interest to us, an SU(2) analysis is 



not enough. We thus study three body decays in the framework of the full SU (3) group. 

The analysis of three body final states involves several subtleties and technical compli- 
cations. We have developed methods to overcome these difficulties that are of more general 
applicability than just the KsK$Ks mode. As concerns our final results, we find that cur- 
rent experimental data give no constraint on the CP asymmetry in B — > K$KsKs using 
only SU(3). It is possible, however, that future measurements of branching ratios of a few 
additional three body modes, together with an improvement in the constraints on a few 
other, will lead to useful constraints. (When we add a rather mild dynamical assumption 
to our SU(3) analysis, we do obtain a bound with present data. The experimental range of 
the CP asymmetries is consistent with this bound.) 

The plan of this paper is as follows. In section |H] we introduce formalism and notations 
that are specifically suitable for three body decays. In section ITTTI we explain the principles of 
how to obtain SU(3) relations that constrain the CP asymmetries in three body decays. In 
section llVl we focus on the mode of interest, B — > KsK$Ks, and give a few concrete examples 
of amplitude relations, as well as a Table that allows one to derive all relevant relations. 
We conclude in section |V| Technical details are further discussed in two appendices. In 
appendix El we derive the relations between the CP asymmetries and the parameters that 
we define in section El and we justify the approximations that we use. In appendix iBl we 
describe the techniques that we developed to deal with the complicated SU (3) decomposition 
of the decay amplitudes. Appendix contains a list of relevant branching ratios. 
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II. NOTATIONS AND FORMALISM 

In this section we show how to modify and generalize the analysis of ref. (^| so that it 
can be applied to three-body decays. 

Unlike two body decays, the final state in three body decays is not uniquely determined by 
the identity of the final mesons. Additional quantum numbers (for example, the momenta) 
are needed to specify the state. We use abstract vector notation, e.g. A KsKsKs , where the 
vector index runs over all possible values for the quantum numbers, to describe the various 
states. The total decay rate is given by 

T(B° - K S K S K S ) = \\A KsKsKs f . (3) 

This equation defines the normalization of the decay amplitudes Af. The norm in the right 
hand side of eq. (J3J) represents a sum over all possible final states. If we choose to describe 
the different final states using definite linear momenta, the norm is actually calculated by an 
integral over all momentum configurations. We stress that the norm is the same, no matter 
which basis we choose to span the final states with. 

In order to derive SU(3) relations, we choose to span the final states in a basis with 
definite linear momenta. Our convention is that the order in which we write the three final 
mesons corresponds to their momentum configuration: 

\MiMjMk) = |M i (p 1 )M 7 (p 2 )M jk (^ 3 )) • (4) 
We further define symmetrized states, \S(f)), as follows: 



|S(MiMiMi)) = \MxMxMt) , 

\S(M 1 M 1 M 2 )) = —/= (|MiMiM 2 ) + \M X M 2 M X ) + \M 2 M X M X )) 
V3 

\S(M 1 M 2 M 3 )) = -±= (|MiM 2 M 3 ) + \M 2 M 3 M X ) + \M 3 M X M 2 ) 
V6 

+ |M 3 M 2 Mi) + \M 2 M X M 3 ) + \M X M 3 M 2 )) 



(5) 



In eq. (J5J), M x , M 2 and M 3 stand for different mesons. 

Focussing on the mode of interest to us, namely -B decay into a final \KsKsKs) state, 
we note it can proceed via any of the three As = — 1 transitions whereby a B° meson decays 
into \K°K°K°), \K°K°K°) or \K°K°K°). Owing to the symmetry of the \K S K S K S ) state 
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under exchange of any two of the final mesons, it can only come from the totally symmetric 
combination of the three states, [^(i^ ^ ^ )). Neglecting CP violation in the neutral kaon 
mixing (the experimental measurement of Ek guarantees that this approximation is good to 
£>(1(T 3 )), we have _ 

(K s K s Ks\S(K°K°K°)) = x/l^r , (6) 

V 8 \ V cs Vcd\ 

for every set of values for the momenta Pi,P2,P3- There are two additional combinations of 
\K°K°K°), \K°K°K°) and \K°K°K°) which are orthogonal to \S(K°K°K )). However, 
since the projection of these combinations on \KsKsKs) is zero, we can write 



Ak s k s k s = V3/8 [(V; s V cd ) I \V: s V cd \\ A s(K0K01 t 0) . (7) 

Within the Standard Model, the violation of CP is encoded in the complex phases of the 
CKM elements. It is therefore convenient, for the purpose of discussing CP asymmetries, 
to have the CKM dependence explicit. Following the discussion above, we thus write the 
B° — > KsK$Ks decay amplitudes as follows: 

■Ak s K s K s — ( VcbVcs a s(K°K°T< ) 

+ v* b v us ^o^o)) Vsjs [(v; s v cd ) i \v; s v cd \] . (8) 

Here, and for all other processes discussed below, the amplitudes for the CP-conjugate 
processes, B° — > f, have the CKM factors complex-conjugated, while the a^'° factors remain 
the same. 

Generalizing j^, we introduce a parameter £: 

p _ [K^Kal a S(K°K°K°) ' a s(K"KQK°) , , 

= \V*V I \\o c II 2 ' ' ' 

where the dot product of complex vectors is defined by X ■ Y = ~Y^ V X* Y v . Another useful 
parameter, is defined as follows: 

\l\ = ■ do) 

\ v cb y cs\ \\ a s(K°K°K0y\ 

We have 

|£| \®S{K°K°K°) ' ®S(K°K°K°)\ < ^ , , 

171 ~~ \\n c II • llfl" II ~~ 

I? I 11 S(K°K°K°)" 11 S(K°K°K )" 

The parameter |£| is the one which can be constrained by SU(3) relations, and that would 
lead, through eq. (fTTj) . to a constraint on 



The case of two body decays [l| constitutes a specific example of our more general nota- 
tion ©, where the vectors are simply one dimensional and, as can be seen from eq. ([lip. 
|£| = |£|. The way in which £ of ref. [2| is defined differs, however, by a weak phase factor: 
f (ref. P)= e*rf(eq. E)- 

Before concluding this section, we introduce one more definition. Experiments often 
measure charge-averaged rates, 

r(B - /) = \ [r(s° - /) + r(5° - 7)] , (12) 



where / is the CP-conjugate state of /. For CP eigenstates, / = ±/. When a single weak 
phase dominates, the CP-conjugate rates are equal, T(B° —>•/) = r(5° — > /), and there is 
no reason to make a distinction between T(B — > /) and T(B° — > /). 

III. CONSTRAINING THE CP ASYMMETRIES 

As we show in appendix 1X1 we can write, to first order in 7£e(£) and Zm(£), 

-Sk s k s k s ~ S^Ks = 2 cos 2/? sin7 7^e(£) , (13) 
Ck s k s k s = -2 shi7Zm(£) . (14) 

The significance of the parameter £ is that it encodes all hadronic physics that affects the 
deviation of —Sk s k s k s from sin 2(3 and of Ck s k s k s from zero. The other parameters, (3 and 
7, are weak phases that can be determined rather accurately from other measurements. We 
learn that if we are able to put an upper bound, |£| < |£| max , we will obtain an unambiguous 
test of the Standard Model CP violation by asking whether the relation [^| 

[{Sksksks + W/cos2/3f + C 2 KsKsKs < 4sin 2 7(|£r x ) 2 (15) 

is fulfilled. As mentioned above, the parameter that appears in the SU(3) relations is 
actually In this and the next section, we assume that SU(3) is exact, and use it to 
constrain |£|. 

In order to constrain |£| we consider As = decay amplitudes and write, using our vector 
notation, 

A f = V: b V cd ¥ f + V: b V ud b u f . (16) 
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SU (3) relations lead to amplitude relations of the form 



->q 

a S(K°K°K0) 



u or C 



(17) 



/ 

Most generally, in the sum over /, states which are permutations of each other are treated 
as different states (for example, K + K~tt° and K + tt°K~ are different states). However, since 
the state \S(K K°K )^ is completely symmetric, the strongest constraint is obtained from 
summing only over completely symmetric terms. Making this choice for eq. (JT7j) allows us 
to rewrite it as follows: 

^(KOKORO) =J2 X f %(/) • ( 18 ) 
/ 

The X/s of eq. (fTHj) are related to the X'/s of eq. (fTTj) by symmetry factors. Taking the 
norm of eq. (|17|) needs to be done with care: the sum can involve states with different 
symmetry properties, and the corresponding norms have different meanings. On the other 
hand, there is no ambiguity in taking the norm of eq. (fTSj) . Consequently, we can write 



/ 



°s(f) 



(19) 



We provisionally assume, for simplicity, that As = ±1 decays are dominated by the 
a c terms (this assumption is justified if |£| is small, see eq. (fTU|) ). while As = decays 
are dominated by the b u terms. (Below we obtain our constraints without making these 
assumptions, in a fashion similar to |9j.) Then the amplitudes are related to the decay rates 
by 



\v; b v cs \ 



"s(k°kok°) ~ V / (8/3)r(5 K s K s Ks) , (20) 
\V: b V ud \\b u s(f) \< ^T{B ^ f) . (21) 

The inequality in eq. (}2T| comes from the fact that we consider the symmetrized state, 
rather then a generic state, for which |V^ 6 \4d| 
(J20j) and ljZI). we get 



T(B -»• /). Combining eqs. (JEU), 



< 



v.. 



V, 



ad 




8/3) T(S - K S K S K S ) 



(22) 



We now proceed_without making the assumptions of cfc- and ^-dominance which led to 
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eqs. pnj) and (jZIj) 9]. Instead of £ , we constrain a new parameter, 



defined by 



I I I I 

Vus ||^c6^cdag {x o K 07fO)+KTi)^ / "<i^( ir Oif07f())|| + \\ VcbV *d a s(K0K0K0)+ V ub V ud^( K K 0K0)\ 



ud 



,+v: h v us di, 



(23) 
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The numerator and denominator of 



are related to charge-averaged rates: 



V cb V c d a S(K°K°K°) VubVud a s(K°K°K°) 



V* V n c 4- V* V TP 1 

v cb y cs U S(K°K K°) ~ v ub v us^ s ^ K o K o K 0) 



<2^|x / |Vr(s-./)j , 

I T/ T/* :7c -4-1/ I/* 

"r \\ VcbV cs u S(K°K t) K ) ' VubV us U S(K K°K ) 

= (16/3)r(S -> K S K S K S ) . 

2 



(24) 



(25) 



Using the measured charge-averaged rates, a constraint on 
further assumptions. 



is obtained without any 



The 



and £ parameters are related as follows: 



VusVcd 
VcsVud 



+ f + 2cos7^e 



VcsV ud 



1 — 1 2 



(26) 



1 + |f I + 2cos7?ee(£) 

The relation (J26)) is a generalization of the relation in j^, eq. 14]. 1 It has the important 



property that for A 2 < 
allowed range, |f| < If 



< 1 we get a constraint on f , for any £ (of course, within the 
see eq. flUJ)). Since we do not know the value of f, we should 
consider the weakest constraint, which corresponds to 7£e(f) = |f| (the (V us V c d) / (V cs V u d) 
term is experimentally known to be real to a good approximation). We show in Fig. ^the 



relation between the upper bound on 



and the resulting upper bound on f , for three 



ranges of the weak phase 7. The weakest bound, which corresponds to 7£e(f) = <M and 



7 = 0, is the curve 



f| — A 2 )/(l + Ifl). Note that the translation from 



non-linear, a point which was not stressed in jjj, although it is true there as well. 



to f is 



IV. 5i7(3) RELATIONS FOR K S K S K S 

The simplest way to find 577(3) relations is to express the decay amplitudes using in- 
variant SU(3) reduced matrix elements. While the number of SU(3) independent reduced 
matrix elements in three body decays is quite large, a significant simplification is obtained 



1 Ref. Q uses the rates T(B a — > /) to define £ , while we use the charge-averaged rates T(B —*■/). If we 
leave 7 unconstrained, the resulting upper bound on |£| is the same, but for 7^0, our expression gives 
a stronger bound, as can be seen in Fig. 2] 
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FIG. 1: The upper bound that can be placed on |£| as a function of the upper bound on |^|, 
according to eq. (|26[). The three curves correspond to different ways of treating the weak phase 
7: 7 = 59.0°, the experimental central value (solid black); 7 £ [35.9°, 80.1°], the 3a range Q| 
(dark-gray dashed); 7 unconstrained (light-gray dotted). 

by the fact that we only consider completely symmetric final states. In particular, for generic 
/'s we have 40 independent reduced matrix elements, while for <S(/)'s there are only 7. 

By scanning over all possible contractions of the relevant SU (3) tensors, we are able to 
obtain all matrix element relations in a systematic way, avoiding the need to discuss SU(3) 
properties of tensor products. We give more details on the calculation in appendix [BJ 

The main results of our work are summarized in table where we list the dependence of 
the symmetrized three body decay amplitudes on SU (3) reduced matrix elements. We give 
here only B° and B + decays, but it is straightforward to add B s decays in a similar way. 
We stress that, since this table includes only totally symmetric states, it is only applicable 
to constrain the totally symmetric final states such as KgKsKs- 

A simple examination of tableQreveals that all the relations of the form of eq. (|18|) involve 
at least one As = decay into a final state with an rj$ meson. In the exact SU(3) limit, this 
corresponds to a state with a final rj meson. We would like to emphasize two points in this 
regard: 

1. The use of SU(3) relations involving amplitudes with final r] 8 and/or r\\ mesons was 
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S(f) 


Ay 


A 2 


A 3 


A 4 


A 5 


A 6 


A 7 


S(K°K°K a ) 


1 




















S(K+K-n°) 





1 

















S(K K tt ) 








1 














S(tt+tt-tt°) 











1 











S{K+tt~K°) 














1 








S(K'tt+K°) 

















1 





S{K+K-tt+) 




















1 


S(K+K°ir°) 





-1 





1 


1 

72 








S(K°K°tt+) 








-V2 


V2 





-1 


1 


S(K°K° m ) 
S(K+K~ m ) 

5(7T+7T — T]g) 

S(K+K°r]8) 

S (lT + IT T]g) 


1 

73 

1 

73 

1 

73 







1 

73 



1 

73 


l 

73 
2 

73 
2 

73 



y§ 




2 
73 

1 

73 

1 

73 







/§ 
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76 

1 

73 






1 

73 









5(7r+7r+7r — ) 




















V2 


5(7r+7r°7r°) 




















1 

72 


S(n n° V8 ) 
S{n° r/gr)g) 
S(ir+ri$ris) 


1 

76 





ATT 

272 
3 

-2 
3 




272 
3 

-2 
3 


i 

76 
5 

372 

4 

3 





72 
3 





72 
3 





1 

72 


5(7r 7r°7r°) 











/I 











s{mmm) 


1 

72 


72 
3 


272 
3 


1 

72 












TABLE I: SU(3) decomposition of As(/)- The different blocks refer to different degrees of sym- 
metrization needed for each state. 

recently criticized in ref. on the basis that SU (3) breaking effects in this system are 
large. The phenomenological value of the octet-singlet mixing angle is sin# w 0.27 
SU(3) breaking in the decay constants is parametrized by 2(/ s — f q )/{f s + fq) ~ 0.22 



|17l |. The breaking effects are thus consistent with our estimated accuracy of (9(0.3), 
and are not 0(1), as suggested in ref. jlti| . 

2. If the relevant decays into final states involving 7/ are found to be enhanced compared 
to the corresponding states involving rj, then the effects of the octet-singlet mixing on 
our results may be significant. These effects can, however, be taken into account by 
using both 77 and 7/ data, in a way similar to ref. jjj. 

We now present several interesting specific relations. Note that, since we do not know 
the values of the strong phases, in deriving our bounds we must add the various As = 
amplitudes constructively, see eq. (|19|). This conservative procedure may weaken the bound 
considerably. Therefore, relations involving a smaller number of As = amplitudes are 
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more likely to give strong bounds. 



A. A Single As = Amplitude 



There is no amplitude relation of the form of eq. (fTHj) that involves only a single As = 
decay amplitude of B° or B + . Such relations would have had the potential to lead to a 
tight constraint. Note that, in general, we would still get an upper bound on, rather than an 



estimate of, 



The reason is that we consider the symmetrized final state while experiments 
measure non-symmetrized final states. However, there exists such a relation involving B s 
decays: 



a B d ^S(K°K K°) 



This relation is, in fact, due to the U-spin subgroup of SU(3) and it holds for the non- 
symmetrized states as well. Since in B s decays the \KsKsKs) state can only come from 
\S(K°K°K )^, in a way similar to eq. (J7J), eq. (|27|) implies relations between the B d — > 
KsK$Ks and B s — > K$KgKs decay amplitudes, V cs V* d a?B d -*KsK s Ks = 

v: s v cd 

b B s ^K s KsK s ' 

leading to 

V us 



V, 



ltd 



K S K S K S ) 



T(B S 

T(B d -> K S K S K S 



(28) 



B. Two As = Amplitudes 



We find a single amplitude relation involving only two As = amplitudes: 



(29) 



The fact that we are interested only in symmetrized states is helpful here in yet another 
way. Let us write 

1 



(S(K K°X)\S(K S K S X)) = -(S(K°K°X)\S(K L K L X)) 



V2' 



(30) 



(X here can be any meson except K°, K°, K$ or K^.) Since in B decays the \S{KgKsX)) 
and \S(KlKlX)) states can only come from an \S(K°K°X)} state, we can write, similarly 
to eq. ©, 



A 



S(K S K S X) 



-A. 



S(K L K L X) 



V2 



A 



S(K0R°X) 



(31) 
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Consequently, the relation ()29|) leads to the following relation, which is more practical from 
the experimental point of view: 

^S(K«K«K») = ^ b S(K s K sm ) ~ ^ b S(K s K s ^y ( 32 ) 

(K$Ks can be replaced by K L K L .) There is yet no measurement of the modes in eq. (|3~2|). 
C. Three As = Amplitudes 

We find several amplitude relations which involve three As = amplitudes, for example, 

At present, the branching ratio B(7r 7r°r]) is not yet constrained, while B(K + K~7t°) and 
B(ti + 7t~7t°) have rather weak upper bounds. For the relation (}3*3*j) to become useful, the 
branching ratio of the first mode must be constrained, and the bounds on the latter two 
must be improved. 



D. Measured As = Amplitudes 

There are relations which involve only modes which have been measured. Given the 
experimental data in appendix El the strongest bound is obtained by using the following 
relation: 



_ hi - \Flh q + \/2 h q 



(34) 



Using the definition 




<0.22W-| W"/ '+,/ Wn+ir-y) / B(ttW 



B(K S K S K S ) V ®(K S K S K S 



(35) 



B{K S K S K S ) V B(K S K S K S) 
As explained above, we take here the SU(3) limit in replacing r] 8 with 77. We also use 
B(S(K°K°n + )) =2B{S{K s K s it+)). 



is too weak to bound |£| and 



We see that the strongest bound we can currently put on 
the CP asymmetries. It is possible, however, that an improvement in experimental data, as 
well as measurements of additional modes, will eventually lead to a significant bound. 
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E. Dynamical Assumptions 



One can use simplifying dynamical assumptions and neglect the effect of small contribu- 
tions from exchange, annihilation, and penguin annihilation diagrams Q|. Practically, this 
means that all reduced matrix elements in which the spectator (the B triplet) is contracted 
with a Hamiltonian operator are put to zero. More details are given in Appendix [Bj 

Such a simplification does lead to new relations. Most notably, there is now an amplitude 
relation involving a single As = mode: 



®%K0K0J?0\ ~ ^^L^q-^o + v • (36) 



3 /B(A>AV+) £ „ 2Q ^ ™ < „ 31 (37) 



This relation leads to the following upper bound 

7 < n 2_ 

" ' "V 2 \jB(K s K s K s ) 

The observed CP asymmetries are well within this bound. We conclude that to uncover a 

signal of new physics with our methods will require improved experimental data. 



V. CONCLUSIONS 



We use the approximate SU(3) flavour symmetry to constrain the SM pollution and 
the CP asymmetries in three body decays. This is an extension of previous works that 
considered two body final states. One important difference is that two body final states are 
entirely defined by the identity of the final mesons, while in three body decays additional 
quantum numbers (such as momenta or angular momenta) are required to characterize the 
final state. In the absence of an experimental spatial analysis, the measured quantities are 
always a sum over all possible final states. On the other hand, the SU(3) relations hold 
for each final state separately. The application of £77(3) relations to three body final states 
should therefore be done with care. 

The case of B — > KsKgKs decay is special since the final state is symmetric under the 
exchange of any two mesons. We showed how this leads to a significant simplification in 
the SU(3) analysis, allowing us to consider only final states with the same symmetry and 
considerably reducing the number of independent SU (3) reduced matrix elements. 

Still, decomposing the decay amplitudes for three body final states into reduced matrix 
elements with well defined ££7(3) transformation properties is a difficult task in terms of 
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group theoretical calculation. Since, however, we are eventually interested only in relations 
between physical decay amplitudes, we were able to scan systematically over all possible 
reduced matrix elements and find an independent subset of them. Using this method, our 
reduced matrix elements bear no clean SU (3) interpretation, but their guaranteed indepen- 
dence is all that matters for the task of finding amplitude relations. The same method can 
be used to simplify other SU(3) calculations where the goal is obtaining physical amplitude 
relations. 

Whether a numerical upper bound is achieved (and whether this bound is strong enough 
to be in conflict with a measured CP asymmetry) depends on the available experimental 
data. Currently, no such bound can be obtained with no additional assumptions, and the 
bound which is obtained when additional dynamical assumptions are used is not strong 
enough to be in conflict with the measured CP asymmetry. However, our work shows which 
new measurements have the potential to lead to a constraint. 

The hope is that, given more and better experimental data, three body decays and 577(3) 
relations will provide us with an additional unambiguous test of the SM mechanism of CP 
violation. 
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APPENDIX A: THE f DEPENDENCE OF S KsKsKs AND C KsKsKs 

In order to derive the dependence of CP asymmetries in B — > K S K S K S on £, we work 
in the basis of states with definite angular momentum between two of the three final Kg 
mesons. While it is difficult to write SU (3) relations in this basis, the advantage of using 
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it is that any basis-state is manifestly CP even (In other bases, final states which are 
related by CP would have related decay amplitudes such that only the combinations which 
correspond to CP even final states have a non-zero amplitude. The discussion, however, is 
much simpler in the basis we choose.) 

We denote the various components of the vector A KsKgKs in this basis by A u where v 
is some collective index which runs over all possible final states. We define for every final 
state the parameter 

K = e** ^ > ( Al ) 
where 4>b is the phase of the B° — B° mixing amplitude (see the review on CP violation in 
jlfll ] and, in particular, eq. (58)). 

Currently, the experimental time dependent CP asymmetry is measured with no distinc- 
tion between various final K$KsKs states. The expression for the measured CP asymmetry 
involves therefore a sum over all final KsK$Ks states: 

AksKsKs {t) = ErM-.i + EW^,) • ( A2) 

V V 

Writing 

•A Ks K s K s {t) = S KsKs k s sin(Am t) - C KsKsKs cos(Am t) , (A3) 
we use the definition flAlj) to get 

£|4,| 2 2ImA„ 

SKsKsKs = £IAf (l + IA.I 2 ) ' (A4) 
EI^| 2 (1-|A.| 2 ) 

° KsKsKs = EI^I 2 (i + |A,l 2 ) ' (A5) 

V 

We can make the CKM dependence of each amplitude explicit: 

K = (V; b V cs al + V: h V us at) [(Ok) / lO^ll = K (1 + iu e 17 ) , (A6) 

where 7 is the weak phase between V* b V cs and V* b V us , and ^ v therefore contains only strong 
phases and is defined by 
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V ub V us 



v*v 



(A7) 
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Assuming that |^| is small for every v (we justify this assumption below), we expand (|A4|) 
and (jA5|) to first order in \^ v \: 

El^fV j:\AH 2 Vu Tie (^) 

Sk s k s k s = sin 2(3 v — + 2 cos 2(3 sin 7 2 , ( A8) 

1^ IA/1 \ A v\ 



E\K\ 2 lm(£u) 

Ck s k s k s = -2 sin 7 - 2 • (A9) 

E \ A t\ 

V 

At this point, the fact that all final states are CP even plays an important role as it dictates 
that r) v — 1 for all v. Switching now to vector notation we have 

Y^\KX ! = ITO 2 Wt , (Aio) 

V 

£ l^l 2 Ke &0 = IOk O^l ^(3° • 3") , (All) 

^|^| 2 Jm(^) = I^KsK^KslXm^-o") . (A12) 

Using the definition Q we therefore get eqs. ()13|) and (|14|) to first order in £. 

We still need to justify the expansions (jA8|) and (jA9|) . In making it, we assumed that for 
every v we have < 1. Unlike in two body decay, the mere smallness of |£| is not enough 
to validate this assumption for every v. Nevertheless, we show next that the smallness of 
|£| does guarantee that the branching ratio of final states in which \^ v \ > 1 is constrained by 
). The terms that are omitted in writing eqs. (jA8|) and (jA9j) are therefore of C?(|£| ). 

Our starting point is the definition of eq. (fTUjl . leading to 

o \V* V I 2 ll«"ll 2 \V*V I 2 II(t"II 2 ■* 

We divide the index i/ into two groups: The group S in which are small and the 
expansions (|A8|) and ()A9|) are justified, and the group S in which |^| > 1 and the expansions 
are not justified. We are interested in the sum over S only. We have 

E_\v: b v us \ 2 K\ 2 TAv: b v us \ 2 K\ 2 + £ I^K/KI 2 

2< — — 2 = |e| 2 + 0(|e| 3 ). (A14) 

However, in the group S, where (V^V^a^l < |V^* b \4 s a"|, the full amplitude \A U \ 2 can be at 
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most 4 |V^* 6 Ks «"| 2 - We therefore find that 



S T< — 2 < 4|e| 2 + 0(|e| 3 ) • (A15) 

V 

One can easily see that omitting these small l^gs] 2 terms from the sums in the expres- 
sions (|A4J) and (JA5j) corresponds to omitting terms of 0(\£\ ) in the expansions ()A8|) 
and (JSnj). 



APPENDIX B: THE SU(S) ANALYSIS 

Finding SU(3) amplitude relations can be done systematically using tensor methods. We 
write down the (tt, K, rj) meson octet as 



V 



7T 

K~ 



7T+ K + 

K° 



\ 



(Bl) 



and the B meson triplet as 



We combine the As 



18. 



20j: 



(B 3 )< = (B + B d b) ■ (B2) 
and As = — 1 Hamiltonian operators into three rank 3 tensors 



my) 



k 



^0 0^ 



VV 












If 



W 




3A^ 
3A^ 



'Xt 



^ ^a; o 







A^ 



A* 



A 



-K 






Uxi 









/ 



A* 



-AS 








-2X d q 










\ / 



V 



3A* 
-A* -AJ 







-2X s q 



II 



(B3) 



(B4) 



(B5) 



where A^ 



v; b v gq , 
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Usually, one proceeds by combining all mesons into irreducible representations of SU(3) 
and contracting the Hamiltonian operators in all possible ways. This would require a large 
amount of multiplications of irreducible representations. Instead, we obtain a set of inde- 
pendent reduced matrix elements by summing systematically over all possible permutations: 



Permutations p \ i\ , 12, 13, 14, is=l 



+4 (^t^t^t^U^tf (B6) 
+A P 15 (Ps)U P ^ P ^U B ^ H ^ 

We remind the reader that the order of the final states mesons is important as it corresponds 
to different momentum configurations. 

All together there are 120 permutations and therefore 360 parameters (A3, A~, A^g). 
However, these 360 free parameters appear in only a small number of combinations which 
correspond to independent reduced matrix elements. Automating the calculation, it is 
straightforward to obtain the set of such independent combinations. Using Young diagrams, 
we verify that we obtain the correct number of independent reduced matrix elements. We 
get 40 independent reduced matrix elements: 10 for H3 in the Hamiltonian (in other words, 
there are 10 reduced matrix elements A3 where ri stands for the 10 different representations 
in 8 x 8 x 8 which have non-zero reduced matrix elements involving H 3 ), 12 for ifg, and 18 
for if 15. The predictive power of the SU(3) symmetry is manifest when one realizes that 
there are 95 As = decays of B + and B° decays to which we can relate the 3 As = — 1 
modes of interest (B° -> K°K°K°, K°K°K° and K°K°K°). 

In this work, only totally symmetric final states play a role. This situation can be used 
to simplify the analysis. We replace the combination (P 8 )^ i (P 8 )^ 2 (P 8 )^ 3 in eq. (|B6|) with 
the symmetrized combination: 

(-fs)^^)^^^)^ + (PsT^iP^Zz^P^Zi 

Then the number of independent matrix elements is reduced to 7. The predictive power is 
maintained, since there are 20 symmetrized As = modes to which we relate the single 

18 



B° — > S(K°K°K°) decay amplitude. Note, however, that the symmetric states are not 
properly normalized and so their normalization needs to be introduced by hand (see eq. (JHj)). 

Table H] lists the independent reduced matrix elements. The names of the matrix elements 
(Ai, A2 etc.) bear no significance. The horizontal lines divide the table according to the 
form of the symmetrized final state. The symmetrized states in the table are all normalized. 

In this work we also consider a simplifying dynamical assumption by which we neglect 
the effect of exchange, annihilation and penguin annihilation diagram la] . The implemen- 
tation of this assumption is straightforward in our calculation. One should just drop all 
permutations in which the 5-meson triplet is contracted with the Hamiltonian operator. In 
other words, one takes eq. (jB6|) and drops from the sum all permutations in which p4 = 4 
or p4 = 5. When this procedure is applied to the symmetrized combination (|B7|) . there are 
only 5 independent reduced matrix elements. 



APPENDIX C: EXPERIMENTAL DATA 



We quote experimental data relevant to three pseudoscalar final states. Measurements 
where resonant contributions are removed from the sample are denoted by (NR). The cur- 
rently measured As = ±1 modes are jjjj]: 



B(K S K S K S 
B(K + it + it- 
B^ nr \K + u + it- 
B(K + K-K + 
B(K + K S K S 
B(t]K + ti- 
B{K\ + 7i- 
B(K + 7i-7i° 
B(K + K-K° 



(5.8 ± 1.0) x 1(T 6 , 
(53.5 ±3.5) x 1(T 6 , 
(4.9 ± 1.5) x 1(T 6 , 
(30.1 ± 1.9) x 10~ 6 , 
(11.5 ± 1.3) x 10" 6 , 
(33.4t«) x 10- 6 , 
(44.9 ±4.0) x 10" 6 , 
(35.6j**) x 10- 6 , 
(24.7 ± 2.3) x 10~ 6 . 



(CI) 
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The currently measured or constrained As = modes are 



QSQ: 



B(K + K~ 



7T 



B(K s K s n + 
B(K + K°n° 
B{K°K-7i + 
B{K + K-7i° 
B{K + K°7i- 



IX TX 7T 



= (16.2 ±2.5) x KT 6 , 
= (16.6ltS) x KT 6 , 

< 6.3 x 10~ 6 , 

< 3.2 x 1(T 6 , 

< 24 x KT 6 , 

< 21.0 x 10~ 6 , 

< 19 x 10" 6 , 

< 18 x 10" 6 , 

< 7.3 x 10~ 6 . 
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